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SUMMARY 


This  report  describes  a procedure  for  identifying  ice  jam  locations  based  on  hydraulic  con- 
ditions in  the  Connecticut  River.  The  hydraulic  analysis  consisted  of  establishing  water  sur- 
face profiles  using  Computer  Program  HEC-2,  developed  by  the  Hydrologic  Engineering  Cen- 
ter, Corps  of  Engineers,  Davis,  California.  Continuity  of  surface  ice  discharge  was  assumed 
in  the  river  reach  studies,  and  this  allowed  computation  of  the  ice  concentration  at  various 
sections  to  determine  the  clogging  potential  of  the  ice  floes. 

Once  conditions  for  clogging  are  identified  at  critical  sections,  a Froude  number  criterion 
is  evaluated  to  determine  ice  block  stability,  and  the  type  of  ice  jam  can  be  identified.  When 
clogging  of  ice  floes  at  a particular  section  is  identified,  depending  on  the  value  of  the  critical 
Froude  number  F and  the  flow  Froude  number  F0I  a jam  forms  by  juxtaposition  of  the  ice 
j/rcr  > F0 ) or  by  underturning  and  stacking  of  the  ice  (Fcr  < F0). 

The  number  of  potential  ice  jam  locations  for  fragmented  ice  floes  decreases  as  the  river 
flow  increases.  This  identification  is  quite  important,  as.  it  implies  that  there  are  more  po- 
tential ice  jam  sites  at  low  flows  because  of  surface  and  subsurface  constrictions  that  arc 
quite  evident  at  low  water  levels.  As  the  water  level  rises,  these  obstructions  are  reduced, 
allowing  greater  surface  conveyance  of  ice. 

The  designs  of  manmade  channels  and  river  modifications  are  changing  to  reflect  the  con- 
ditions in  which  the  natural  channels  are  formed.  Channels  with  manmade  bends  are  being 
designed  to  accommodate  fish  habitats  as  well  as  to  reduce  potentially  high  sediment  trans- 
port. Reaches  with  trapezoidal  sections  to  replace  oxbow  sections  are  being  discouraged 
and  conformity  with  more  existing  natural  conditions,  such  as  keeping  the  oxbows  but  im- 
proving the  conveyance  capacity  and  stabilizing  the  banks,  is  being  stressed.  For  design  of 
such  channels  in  cold  regions,  an  ice  jam  location  check  such  as  the  one  described  in  the 
report  would  ensure  that  unfavorable  conditions  conducive  to  ice  jamming  are  avoided. 
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ANALYSIS  OF  POTENTIAL  ICE  JAM  SITES  ON  THE 
CONNECTICUT  RIVER,  AT  WINDSOR,  VERMONT 


D.J.  Calkins,  M.S.  Hutton  and  T.L.  Marlar 


INTRODUCTION 

The  location  of  ice  jams  along  a river  system  and  the  prediction  of  their  occurrences  have  previ- 
ously been  studied  through  review  of  historical  records.  The  systematic  analysis  of  the  hydraulic 
conditions  and  river  ice  characteristics  that  initiate  ice  jams  has  been  used  only  rarely.  Often  the 
destructive  end-products  of  ice  jams  have  been  overstressed,  whereas  the  conditions  initiating  them 
have  been  overlooked,  perhaps  because  of  the  spectacular  results  of  massive  ice  piling. 

A joint  project,  therefore,  was  begun  by  the  New  England  Division,  Corps  of  Engineers,  and  the 
U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory  to  study  the  formation  and  breakup 
of  ice  on  the  Connecticut  River  between  White  River  Junction  and  Windsor,  Vermont,  during  the 
1974-1975  season  (Fig.  1).  Past  records  indicated  that  frequent  ice  jams  occurred  in  this  reach  of 
the  river. 

After  initial  review,  it  became  evident  that  the  river  reach  selected  was  too  large  for  complete 
analysis;  therefore,  a shorter  section,  from  the  railroad  bridge  in  Windsor  upstream  to  Sumner  Falls, 
was  chosen  for  investigation.  The  total  length  of  this  shorter  reach  was  approximately  36,000  ft  as 
measured  from  U.S.  Geological  Survey  topographic  maps.  Low-altitude  aerial  photography  was 
used  to  delineate  stream  patterns,  rock  outcrops  and  manmade  structures,  which  are  often  associ- 
ated with  the  initiation  or  formation  of  ice  jams.  Aerial  photography  was  also  used  to  document 
ice  conditions  in  the  study  reach. 


TEST  SITE 

In  August  1974,  an  aerial  photographic  mission  (altitude  2000  ft,  scale  1:3250)  was  accomplished 
over  the  proposed  area,  and  from  the  photography  obtained,  a smaller  region  of  the  river  was  selected 
for  detailed  analysis  (Fig.  1 ).  This  consisted  of  a 6.5-mile  stretch  of  the  Connecticut  River  from  Sum- 
ner Falls  in  North  Hartland,  Vermont,  to  the  southern  tip  of  Chase  Island  in  Windsor,  Vermont. 

A flood  control  reservoir  exists  upriver  from  the  test  site  on  the  Ottauquechee  River  in  North 
Hartland.  The  outflow  from  the  North  Hartland  Dam  is  regulated  only  during  high  flows  or  antici- 
pated extreme  rainfall  events,  so  the  inflow  to  the  Connecticut  River  is  generally  uncontrolled.  The 
Wilder  Dam  on  the  Connecticut  River  in  Lebanon,  New  Hampshire,  releases  water  on  a predeter- 
mined schedule  set  by  the  New  England  Power  Company.  This  release  normally  occurs  every  morn- 
ing, Monday  through  Friday,  to  generate  hydroelectric  power,  and  as  a consequence,  the  water  level 
in  the  test  area  rises  up  to  5 ft  in  late  morning  to  early  afternoon  depending  on  the  discharge. 
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Figure  1.  Location  of  test  site. 


A typical  discharge  hydrograph  obtained  at  White  River  Junction  gaging  station  is  shown  in  Fig- 
ure 2 for  22  October  1974.  The  duration  and  rise  in  the  hydrograph  are  generally  determined  by 
the  anticipated  peak  power  demand  and  water  release  schedules  set  for  the  entire  Connecticut 
River  Basin.  Weekend  hydropower  generation  is  not  usually  required  at  the  Wilder  Dam  or  for  the 
other  upstream  power  dams  unless  unexpected  peak  electrical  demands  are  reached. 

The  other  significant  tributaries  of  the  Connecticut  River  in  the  study  area  are  the  Mascoma  and 
White  Rivers,  entering  at  West  Lebanon,  New  Hampshire,  and  White  River  Junction,  Vermont,  re- 
spectively, with  the  White  River  Basin  being  the  larger  of  the  two.  The  river  gaging  station  for  the 
Connecticut  River  at  White  River  Junction  includes  the  flow  from  the  White  River.  The  Mascoma 
River  is  the  only  river  with  a potentially  significant  ungaged  winter  flow. 

The  flood  control  structure  at  North  Hartland  and  the  many  upstream  mill  ponds  retain  all  the 
ice  produced  in  the  Ottauquechee  River.  The  ice  supply  from  the  Mascoma  River  is  also  limited 
because  of  several  small  control  structures  between  West  Lebanon  and  Lebanon  which  retain  the 
ice.  The  significant  ice  discharge  to  the  study  reach  is  that  generated  in  the  White  River  and  in  the 
Connecticut  River  up  to  Wilder  Dam.  The  Wilder  Dam  retains  all  the  ice  generated  in  the  Connec- 
ticut River  and  its  tributaries  up  to  Barnet,  Vermont. 


AERIAL  PHOTOGRAPHY 

Photography  obtained  for  the  study  of  ice  jams  in  the  Connecticut  River  consisted  of  9-in.  by 
9-in.  vertical  photographs  from  a Zeiss  RMK  15/23  mapping  camera  with  a 6-in.  lens,  70-mm  verti- 
cal photographs  from  a Hasselblad  500EL  camera  with  a 100-mm  planar  lens,  and  70-mm  oblique 
photographs  from  a Hasselblad  500EL  camera  with  80-mm  and  100-mm  lenses.  Films  used  were 
Kodak  Plus-X  aerographic  (type  2402),  Kodak  Tri-X  professional  panchromatic,  and  Kodak 
Ektacolor  professional  film  (type  S).  The  9-in.  by  9-in.  photographs  were  used  to  prepare  a photo- 
mosaic of  the  test  area.  The  70-mm  photographs  were  used  to  document  the  ice  buildup  and  re- 
lated ice  features  throughout  the  winter. 

Stereographic  photography  was  obtained  at  an  altitude  of  2000  ft  above  mean  sea  level  and  this 
resulted  in  a variable  scale  of  photography  because  of  changes  in  terrain  elevation  along  the  river. 

A mean-sea-level  altitude  was  chosen  because  it  provided  a basis  for  scale  correction  and  because  it 
was  not  possible  to  maintain  a constant  distance  above  ground  level  due  to  rapid  terrain  elevation 
changes.  The  scale  corrections  necessary  for  the  river  portion  of  the  photographs  were  minimal 
because  the  elevation  drop  of  the  water  surface  between  the  upstream  section  and  the  downstream 
section  was  less  than  1 5.0  ft. 

To  reduce  the  cost  of  acquiring  aerial  photography,  a small  format  camera  system  was  designed 
by  Thomas  Marlar  and  constructed  at  USA  CRREL.  This  utilized  the  door  of  a Cessna  T-41  air- 
craft and  various  off-the-shelf  camera  components  (Fig.  3).  An  additional  door  was  purchased  and 
modified  so  that  it  could  accept  the  camera  system  and  still  remain  easily  interchangeable  with  the 
standard  aircraft  door.  This  system  is  usable  on  all  Cessna  1-41/172  model  aircraft  and  can  be 
readily  adapted  to  other  high-wing  aircraft.  A few  of  the  ice  Jam  surveillance  flights  were  flown 
using  this  system  to  obtain  vertical  photographs  and  excellent  results  were  achieved. 


CHANNEL  MORPHOLOGY  - INTERPRETATION  OF  AERIAL  PHOTOGRAPHS 

An  uncontrolled  mosaic  of  the  Connecticut  River  from  Sumner  Falls  to  Chase  Island  was  as- 
sembled. The  scale  of  the  working  photographic  prints  was  approximately  1 :3250.  Stream  lea- 
tures  such  as  alluvial  fans,  rock  outcrops,  log  debris,  channel  sinuosity,  bank  characteristics  and 


Figure  3.  Views  of  the  surveillance  cameras  mounted  in  the  door  of  a Cessna  T-41  aircraft. 
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others  could  be  readily  identified  from  low-altitude  black  and  white  photographs.  The  vertical 
photographs  that  are  presented  in  this  report  (Fig.  4)  have  been  reduced  from  these  original  9-in. 
by  9-in.  prints. 

Figures  4a-g  show  photographs  taken  with  the  mapping  camera  during  the  month  of  August  at 
low  river  flow  conditions  (Q  « 1000  ft3 /sec).  The  locations  of  the  cross  sections,  islands,  and 
tributaries  along  the  test  reach  are  presented  in  the  Appendix. 

Areas  in  a stream  channel  where  ice  jams  generally  form  are:  1 ) constrictions;  2)  rock  outcrops 
and  manmade  structures  (bridge  piers);  3)  long,  slow-velocity,  deep  water  pools;  and  4)  shallow  sec- 
tions across  portions  of  the  channel  where  grounding  of  ice  floes  could  be  initiated.  There  were 
three  identifiable  surface  constrictions  in  the  study  reach,  two  occurring  as  a result  of  bridge  piers 
and  the  third  at  Section  1 2.  Section  1 2 occurs  at  a sharp  bend  (~  70°)  which  has  riprap  on  the 
outside  bank  to  protect  the  Boston  and  Maine  Railroad  bed  (Fig.  4d).  All  visible  rock  outcrops 
were  identified  and  they  were  located  between  Sections  7 and  8 (Fig.  4c). 

The  other  potential  ice  jam  sites  occur  between  Sections  18  and  19  and  Sections  5 and  7,  where 
significant  widening  of  the  channel  has  occurred  with  deposition  zones  having  been  formed  into 
sandbars,  visible  at  low  flows  (Fig.  4b  and  c).  The  material  in  these  areas  was  classified  as  a silty  sand, 
while  the  bed  material  found  in  the  main  sections  of  the  river  was  classified  as  a poorly  graded  gravel. 

GROUND  SURVEYS 

During  October  1974,  cross  sections  were  surveyed  in  the  study  reach.  The  sections  were  located 
after  examination  of  black  and  white  aerial  photographs.  Two  additional  cross  sections  were  added 
when  the  actual  surveying  was  being  conducted.  The  survey  was  initialed  at  the  crest  of  Sumner 


b.  Slow-velocity  reach  just  above  Section  IS.  1)  River  has  widened  with  a deposition 
zone;  exposed  bed  material  and  debris  can  be  seen.  2)  Ripple  section  begins  at  the  end 

of  this  long  pool. 


c.  Another  slow-velocity  section  near  Station  6,  opposite  Windsor.  1)  Main  channel 
meander.  2)  Deposition  zones  with  the  bottom  detail  easily  shown  as  light  gray. 


Figure  4 (cont ’d).  Aerial  photographs  of  study  area. 


— til 


a 

9 ‘x  ^ v A 

Wf  JEi?  ® 1 

KjK§| 

m 


■ 


t 


> 

\ 


Falls  with  an  assumed  starting  elevation  of  400.00  ft.  To  convert  these  elevations  to  the  U.S.  Geo- 
logical Survey  (USGS)  bench  marks,  60.6  ft  has  to  be  subtracted.  Temporary  bench  marks  were 
established  along  the  river  beginning  at  Station  1 2.  These  spot  elevations  were  located  one  to  two 
sections  apart,  and  were  established  to  allow  elevation  readings  to  be  taken  on  any  potential  ice 
jam  within  the  study  reach. 

The  actual  surveying  of  the  cross  sections  was  performed  by  using  two  methods:  wading  with 
a surveying  rod,  and  sounding  from  a small  boat.  The  sounding  method  involved  two  techniques, 
depending  on  water  velocity  and  depth.  The  boat  was  positioned  by  being  headed  upstream  into 
the  current.  The  sounding  line  (marked  at  1 -ft  increments)  with  a 2-lb  weight  at  the  end  was 
thrown  slightly  upstream  into  the  current,  and  the  depth  to  the  bottom  was  recorded  when  the 
weight  contacted  the  riverbed.  When  the  velocit>  was  greater  than  3 ft/sec  and/or  the  depth 
greater  than  10  ft,  a heavier  sounding  weight  (~  1 S lb)  was  used.  The  horizontal  position  in  the 
cross  section  was  measured  using  stadia  from  the  surveying  rod.  The  accuracy  at  300  ft  was  ap- 
proximately ±6  ft  because  of  bobbing  of  the  boat.  Proper  alignment  of  the  boat  in  the  cross  sec- 
tion was  maintained  through  the  use  of  the  transit  set  up  on  the  bank  at  one  end  of  the  cross  sec- 
tion. The  distance  between  cross  sections  was  taken  from  the  uncontrolled  mosaic  of  the  black 
and  white  aerial  photographs.  Plots  of  the  cross  sections  arc  shown  in  the  Appendix. 

A continuous  water-level  recorder  was  installed  on  State  of  New  Hampshire  property  in  Cornish, 
approximately  one-half  mile  upstream  from  the  W'indsor  Covered  Bridge,  near  the  boat  landing 
access.  It  was  installed  during  January  (with  assistance  of  personnel  from  the  USGS  in  Concord, 
New  Hampshire,  and  the  Reservoir  Control  Center,  New  England  Division,  Corps  of  Engineers)  to 
record  the  changes  in  water  level  during  the  formation  stage  of  the  anticipated  ice  jam.  However, 
no  ice  jam  occurred  during  the  ice  season,  and  only  peak  levels  due  to  snowmelt  were  recorded. 


HYDRAULIC  ANALYSIS 


The  basic  hydraulic  analysis  was  limited  to  the  determination  of  the  water  surface  profiles  as 
calculated  using  the  Computer  Program  HEC-2,  developed  at  the  Hydrologic  Engineering  Center, 
Corps  of  Engineers,  Davis,  California.  This  method  of  calculation  uses  Bernoulli’s  equation  for  the 
total  energy  at  each  section  and  Manning’s  formula  for  determining  the  head  loss  between  sections. 
Expansion  and  conti action  coefficients  were  assumed  to  be  on  the  order  of  0.3  and  0.1 , respectively. 
The  energy  losses  due  to  the  bridge  piers  were  ignored  because  the  necessary  field  data  for  using 
the  special  bridge  routine  were  not  taken  and  their  effects  were  assumed  to  be  minimal. 

The  premise  for  using  this  program  was  to  establish  the  usefulness  of  the  data  generated  for 
analyzing  potential  ice  jam  locations  in  a river  system.  Also,  this  program  is  available  at  most 
Corps  of  Engineers’  division  and  district  offices.  A contract  with  the  Dufrcsne-Henry  Engineering 
Corporation  in  North  Springfield,  Vermont,  provided  computer  support  in  running  the  HEC-2 
program  on  an  IBM-370  system.  Data  on  the  following  river  characteristics  were  obtained  for  in- 
put to  the  program:  1)  channel  cross  sections,  2)  Manning's  roughness  coefficient  (n)  for  the  chan- 
nel, and  3)  starting  elevations  corresponding  to  three  specific  discharge  values  (i.e.,  a rating  curve 
for  the  downstream  section  of  the  river). 


The  data  for  the  channel  cross  sections  were  taken  from  the  field  survey,  and  the  distances  be- 
tween sections  were  scaled  from  the  photomosaic.  Manning’s  n was  estimated  by  taking  field  ob- 
servations of  the  river  bottom  and  bank  characteristics,  and  comparing  them  with  photographs  of 
similar  channels  with  previously  determined  values  of  Manning's/)  (Chow  1959).  It  is  recognized 
that  many  factors,  such  as  bed  material,  vegetation,  channel  alignment,  discharge,  etc.,  influence 
Manning’s/).  Values  of  n for  the  flood  plain  were  not  considered  because  the  flows  analyzed  for 
this  study  were  contained  within  the  banks. 


9 


0 4 8 12  16  20  24  28  32x1000 

Distonce  Upstream,  ft 

Figure  5.  Water  surface  profiles. 

The  most  difficult  data  to  obtain  were  the  starting  elevations  to  be  used  at  the  downstream  sec- 
tion (Section  2)  corresponding  to  specific  discharges.  The  nearest  continuous  river  gaging  station 
was  located  about  10  miles  upstream.  The  Mascoma  and  the  Ottauqucchee  Rivers  flow  into  the 
Connecticut  River  below  that  gaging  station,  but  above  the  study  reach.  Two  river  flows  of  1 ,000 
and  9,000  ft’/sec  were  gaged  at  the  downstream  section  and  at  selected  points  upstream.  Three 
additional  flows  5,000,  1 7,000  and  25,000  ft3 /sec  were  to  be  analyzed  using  the  HEC-2  program, 
but  there  was  no  stage-discharge  curve  available  for  the  downstream  section.  The  approximate 
elevate  ns  for  these  flows  at  Section  2 were  estimated  based  on  the  stage  discharge  data  from  the 
White  River  Junction  station  and  the  two  measured  flows  at  Section  2,  with  adjustments  in  stage 
height  made  for  the  flows  from  the  Mascoma  and  Ottauqucchee  Rivers. 

A convergence  test  for  the  water  surface  profile  was  performed  lor  a flow  of  9,000  ft3 /sec. 
Cross  section  2 was  selected  as  representative,  and  a distance  of  20,000  ft  downstream  was  used 
as  the  length  of  the  channel.  The  values  of  the  parameters  were:  n = 0.03;  bed  slope  = 0.0028  ft/ 
ft;  and  starting  elevation  = 353.1 . The  computed  water  surface  elevation  at  Section  2 was  357.97. 
The  measured  starting  elevation  was  357.7  for  9,000  ft3  /sec;  since  this  difference  between  the  cal- 
culated and  measured  starting  elevations  is  within  acceptable  tolerances,  the  calculations  proceeded 
using  357.7 

A trial  run  was  performed  using  Manning's  n of  0.040  for  a discharge  of  1,000  ft3 /sec  and 
comparing  the  water-surface  profile  using  an  n value  of  0.030.  The  computed  water  levels  at  Sec- 
tion 1 2 were  357.2  and  357.5  for  n values  of  0.030  and  0.040,  respectively.  This  small  difference 
was  not  considered  significant  for  the  present  project;  therefore,  all  subsequent  profiles  were  cal- 
culated using  n = 0.030. 
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To  assess  possible  errors  in  the  predicted  water  surface  profiles,  the  starting  water  surface  ele- 
vations for  each  discharge  were  run  with  a 1.0-ft  increment  in  elevation  difference.  The  difference 
in  elevation  at  Section  1 2 between  the  two  computed  profiles  for  all  discharges  was  less  than  0.1  5 
ft. 

For  the  low  flow  of  1,000  ft3/scc,  critical  depths  were  calculated  at  Sections  3,  13  and  15.  Fig- 
ure 5 shows  the  water  surface  profiles  for  3 discharges  with  the  starting  elevation  of  358.7  for  the 
flow  of  9,000  ft3/sec.  Flows  exceeding  25,000  ft3 /sec  were  not  considered  because  the  initiation 
stage  of  ice  jams  is  thought  to  occur  at  low  to  moderate  flows.  Detailed  water  surface  profiles  for 
these  discharges  are  on  file  at  USA  CRREL  as  Internal  Report  423  (Calkins  1975b). 


ANALYSIS  OF  ICE  JAM  LOCATIONS 


The  river  sites  where  ice  jams  occur  have  generally  been  identified  through  historical  records. 
Attempts  to  analyze  systematically  the  hydraulic  conditions  and  ice  characteristics  that  initiate 
the  ice  jams  have  been  minimal.  Several  studies  have  concentrated  on  the  equilibrium  thickness  of 
an  ice  jam  once  it  has  formed;  Uzuner  and  Kennedy  (1974)  have  summarized  the  various  methods. 

Nuttall  (1973)  presented  some  basic  work  on  identifying  ice  jam  sites  based  on  surface  ice  dis- 
charge and  a type  of  Froude  criterion.  Calkins  (1975a)  took  the  same  data  used  by  Nuttall  and 
evaluated  the  critical  Froude  number  based  on  block  thickness  and  flow  depth.  This  identification 
scheme  for  the  location  of  jam  sites  was  based  on  fragmented  ice  flowing  downstream  in  a river 
system. 

Continuity  of  surface  ice  discharge  is  assumed  with  no  storage  of  ice  floes  between  sections. 

The  ratio  xc  of  ice  surface  flux  (ft2 /sec)  between  any  downstream  section  and  the  upstream  con- 
trol section  (chosen  for  this  study  as  Section  18)  is 


xc 


w i 


(i) 


where  V = surface  velocity  (ft/sec) 

B = channel  width  (ft) 

C = surface  concentration  of  floating  ice  (0-1.0) 
and  the  subscripts  i and  u represent  any  downstream  section  and  Section  18,  respectively. 

If  one  assumes  the  flow  discharge  is  constant  between  reaches,  = Qu,  eq  1 can  be  rearranged 
to  yield 
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where  A is  the  channel  cross-sectional  area  at  the  respective  sections.  Equation  2 could  have  been 
written  in  terms  of  the  hydraulic  radius  /?h,  as  the  term  A/B  closely  approximates  A/P  where  P is 
the  wetted  perimeter  for  open-water  conditions. 

Frankenstein  and  Assur  (1972)  hypothesized  that  clogging  of  ice  floes  would  begin  when  the 
ice  concentration  reached  0.707,  which  probably  is  a reasonable  estimate.  A jammed  section 
might  be  classified  as  an  advanced  stage  of  clogging,  reaching  a value  of  C = 1 .0. 
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Figure  6.  Critical  Froude  number  based  on  flow  depth  (Ashton  1974a). 

The  leading  (upstream)  edge  stability  of  individual  floating  ice  blocks  was  analyzed  by  a num- 
ber of  authors,  including,  most  recently  Ashton  (1974a,  1 974b)  and  Uzuner  and  Kennedy  (1972). 
Ashton’s  work  is  used  here  where  the  critical  Froude  number  is  calculated  based  on  a combination  of 
moment  equilibrium  and  hydrodynamic  analyses  to  determine  the  values  of  hydraulic  parameters 
for  entrainment  of  the  ice  floes.  The  relation  between  t/H  (ice  thickness  to  flow  depth  ratio)  and 
Fcr  (critical  Froude  number)  is  given  in  Figure  6,  from  which  is  derived  the  relationship 

F = = 7(\-tlH)[tlHf  (3) 

cr  |<?W(1-p7p)1*  [ 5-3(1  -r/w)2]* 

where  V = critical  velocity  for  entrainment  of  an  ice  block  at  the  leading  edge  of  an  accumulation 
g = gravity 
t = ice  floe  thickness 

H = depth  of  flow  at  the  upstream  edge  of  the  cover 
p'lp  = 0.92  for  solid  ice  floes,  neglecting  porosity  effects. 

The  flow  Froude  number  FQ  for  any  section  can  be  computed  if  the  depth  H and  flow  velocity 
V are  known 


The  critical  values  of  the  Froude  number  F can  be  calculated  using  various  ice  block  thick- 
nesses, since  the  backwater  analysis  yields  the  mean  velocity  and  flow  depth  as  well  as  other 
hydraulic  parameters.  It  was  assumed  for  this  analysis  that  the  mean  channel  velocity  is  equal  to 
the  surface  velocity  although  it  is  recognized  that  the  surface  velocity  might  be  in  error  by  5%, 
but  the  assumption  is  within  reasonable  limits. 

The  Froude  number  criterion  described  above  evaluates  the  stability  of  individual  ice  floes 
encountering  the  upstream  edge  of  the  ice  cover.  If  F0  < F , then  the  individual  ice  floes  are 
stable  and  will  not  underturn  and  the  ice  cover  will  lengthen  upstream  by  floe  juxtaposition  until 
the  condition  F0  > Fcr.  When  F0  > Fct  the  ice  floes  will  be  unstable;  underturning  and  stacking 
of  the  floes  will  occur  with  potential  transport  and  accumulation  beneath  the  ice  cover.  The  transport 
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distance  depends  on  local  velocity  conditions,  block  size,  channel  characteristics  and  other  hydraulic 
conditions  (Ashton  1975). 

The  ice  surface  discharge  parameter  xc  indicates  the  potential  reaches  in  a channel  where  the 
congestion  of  ice  floes  is  occurring.  When  the  ice  concentration  equals  or  exceeds  1 .0,  it  may  be 
assumed  that  individual  ice  pieces  have  to  accumulate  by  thickening  or  shoving,  indicating  a poten- 
tial ice  jam  location. 

The  data  used  in  the  following  figures  were  generated  by  the  Computer  Program  HEC-2  except 
where  assumptions  were  made  concerning  ice  thickness.  Figure  7 depicts  the  channel  width  B,  flow 
Froude  number  F0  and  the  ice  surface  discharge  parameter  xc  for  flows  of  1,000,  5,000,  9,000,  and 
25,000  ft3/sec.  For  the  low  flow  of  1 ,000  ft3/sec,  critical  depth  occurs  at  three  sections,  as  indi- 
cated by  F0  = 1 .0.  The  surface  ice  discharge  parameter  xc  is  better  portrayed  in  Figure  8 for  flows 
of  1,000,  9,000  and  25,000  ft3 /sec. 

At  low  flows,  xc  is  highly  variable,  indicating  channel  constrictions  such  as:  Sections  2A  and  4, 
which  are  just  upstream  from  the  bridge  piers;  Section  7A,  where  the  channel  narrows  and  deepens 
just  downstream  from  the  ledge  outcrops;  and  Section  1 2,  which  is  at  a sharp  bend  in  the  river. 

These  features  are  easily  seen  from  the  aerial  photographs  (Fig.  4d-e).  As  the  flows  increase,  the 
parameter  xc  fluctuates  less  and  approaches  the  value  of  1 .0.  This  is  to  be  expected,  since  as  the 
river  stage  rises  within  the  main  channel,  local  sandbars,  minor  obstructions,  and  other  bottom 
features  are  less  dominant,  and  in  this  reach  the  stream  top  width  fluctuates  little. 

To  calculate  the  downstream  ice  concentration,  the  incoming  ice  concentration  at  Section  18 
(upstream  control  section)  is  assumed,  and  then  multiplied  by  xc  to  compute  the  percentage  of 
coverage  by  ice.  For  example,  if  the  ice  concentration  at  Section  18  for  a flow  of  1,000  ft3 /sec  is 
Cu  = 0.60,  then  the  expected  concentration  C-t  at  Section  1 2 would  be  1 .21 , where  xc  = 2.02.  This 
indicates  clogging  of  ice  floes  at  Section  1 2 and  suggests  the  possibility  of  a jam  forming  in  this 
location.  Other  potential  sites  are  indicated  in  Figure  10  by  the  peak  values  xc-  Thus,  for  a flow 
of  25,000  ft3 /sec  at  Section  1 2 (Cu  = 0.60,  xc  = 1 .25),  Cj  = 0.75,  and  jamming  of  ice  at  this  section 
is  not  a problem  because  of  the  increased  available  surface  area  for  ice  floe  transport. 

Figure  9 presents  a frequency  distribution  of  potential  occurrences  along  the  entire  river  where 
the  value  of  C(  has  exceeded  1 .0  for  three  upstream  ice  surface  concentrations  and  three  selected 
flow  rates.  Note  that  the  frequency  of  ice  floe  clogging  occurrences  decreases  as  the  flow  rate  in- 
creases. This  implies  that  ice  jams  from  fragmented  floes  are  initiated  at  low  discharges,  not  as  a 
result  of  high  discharges. 

At  sections  where  Ci  > 1 .0,  the  critical  Froude  number  Fcr  is  next  examined  to  evaluate  block 
stability  at  the  jammed  section.  Figures  lOa-c  depict  the  variation  in  flow  Froude  number  and  the 
critical  Froude  number  along  the  river.  For  flows  of  25,000  and  9,000  ft3 /sec,  F0  is  greater  every- 
where than  F for  assumed  ice  thicknesses  of  1,  2,  and  3 ft.  For  the  flow  of  9,000  ft3/sec,  the  F0 
approaches  Fcr  at  Sections  2A,  4,  7 and  12. 

At  the  low  flow  of  1,000  ft3/sec,  F0  is  less  than  Fcr  at  several  sections;  this  indicates  individual 
ice  block  stability  (no  thickening  of  accumulated  ice)  at  these  sites.  Also,  certain  values  of  Fcr  are 
equal  to  zero  at  Sections  3,6,  11,  13  and  1 5.  This  corresponds  to  grounding  of  ice  floes  on  the 
channel  bottom.  Figure  10a  depicts  this  highly  variable  flow  condition  for  Q = 1,000  ft3 /sec,  as 
indicated  by  the  Froude  number,  implying  ice  block  stability,  instability,  or  grounding  of  the  floes. 

The  potential  ice  jam  locations,  as  indicated  by  the  parameter  xc  for  a discharge  of  1,000  ft3 /sec, 
are  at  Sections  2A,  4,  7A  and  1 2 (see  Fig.  8).  Without  any  knowledge  of  the  upstream  ice  condi- 
tions, one  could  assume  that  potential  ice  jam  sites  are  located  at  the  sections  where  xc  reaches 
maximum  values.  For  an  ice  thickness  of  1 .0  ft,  the  flow  Froude  number  F„  is  less  than  the  critical 
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Channel  parameters,  width  B,  flow  Froude  number  F^  and  ice  surface  discharge  \ CD  = critical  depth. 


Figure  7 (corn’d). 
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Figure  10.  Critical  and  flow  Froude  numbers  along  the  channel  for  three 
ice  thicknesses. 
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Figure  10  (con t’d). 


Froude  number  Fcr  at  all  four  sections  (Fig.  10),  indicating  no  underturning  or  stacking.  If  the  ice 
floes  are  of  sufficient  internal  strength  and  size,  arching  of  the  ice  floes  could  be  expected  at  these 
four  sections.  The  immediate  sections  upstream  of  2A  and  12  have  high  values  of  F0  and  low  values 
of  xc,  indicating  that  clogging  is  minimal  and  that  underturning  of  ice  floes  will  take  place.  Upstream 
from  Sections  4 and  7,  the  ice  will  accumulate  by  juxtaposition  until  the  condition  F0  > Fcr  is 
reached.  Realizing  the  high  variability  of  any  ice  data  and  the  resulting  complicated  flow  patterns 
at  Q * 1 ,000  lt3/sec,  the  F0  values  for  the  entire  reach  between  Sections  4 and  7A  are  below  or  near 
the  critical  values,  and  thus  indicate  a reach  conducive  to  stabilizing  an  ice  cover  of  fragmented 
floe,. 

When  the  river  flow  increases  from  1,000  to  9,000,  then  to  25,000  ft3/sec,  the  difference  be- 
tween the  flow  and  critical  Froude  numbers  increases,  and  the  individual  ice  blocks  become  un- 
stable. The  variability  in  the  ice  discharge  parameter  xc  diminishes  as  the  flow  increases,  because 
of  the  improved  surface  conveyance  of  the  ice,  thus  minimizing  the  potential  clogging  of  ice  floes 
along  the  river. 

Once  the  ice  cover  has  been  formed  at  low  flow  conditions  between  Sections  2A  and  10  and 
during  the  stage  rise  every  day,  the  incoming  broken  ice  from  upstream  will  probably  be  accumu- 
lated beneath  the  cover  because  the  critical  Froude  number  is  exceeded  when  the  flow  is  greater 
than  5,000  ft’/sec  and  the  release  flow  from  Wilder  approaches  10,000  ft3 /sec.  The  constant  ris- 
ing and  falling  water  levels  in  the  reach  every  day  complicate  the  ice  floe  behavior. 


WINTER  FIELD  OBSERVATIONS 

The  formation  of  an  ice  cover  is  complex,  with  shore  ice  often  the  first  to  appear  in  the  quiescent 
zones  along  banks.  Frazil  ice  formation  in  a river  system  is  associated  with  turbulent  flow  and  super- 
cooling of  the  water.  The  temperature  associated  with  nucleation  of  frazil  ice  is  only  a few  hun- 
dredths (~  0.05°)  below  32°F.  Once  the  frazil  crystals  leave  the  so-called  active  zone,  they  tend  to 
concentrate  at  the  water  surface,  and  the  development  of  slush-icc  pans  or  floes  from  the  frazil  ice 
depends  in  part  upon  the  velocity  conditions  in  the  river. 

Frazil  ice  is  developed  in  the  White  River  and  moves  into  the  Connecticut  River  at  White  River 
Junction.  Frazil  is  also  produced  in  the  Connecticut  River,  probably  below  Sumner  Falls.  Frazil 
production  does  not  usually  take  place  in  water  velocities  less  than  2.5-3. 0 ft/sec.  The  average 
stream  velocity  in  the  Connecticut  River  is  greater  than  2.5  ft/sec  for  flows  exceeding  5,000  ft3 /sec. 

The  quantity  of  frazil  ice  production  can  be  large;  for  example,  a basic  heat  transfer  relationship 
has  been  formulated  to  estimate  the  ice  production  in  a river  reach  (Michel  1971 ): 
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where  = discharge  of  ice  (ft3 /day) 

<1/1  = average  heat  extracted  from  stretch  / 
b-  = width  of  open  water  on  stretch  / (ft) 
A/j  = length  of  stretch  (ft) 

L = latent  heat  of  fusion  of  ice  (Btu/lb) 
Pj  = density  of  ice  (slug/ft3 ) 
g = gravitational  constant  (32.2  ft/sec2) 


(Btu/ft2  day) 


If  we  assume  L = 144,  p(  = 1.78,  A/,  = 1.0,  bi  = 1.0  and  To  = 1000  for  an  air  temperature  of  roughly 
20°F,  the  ice  discharge  is  = 0.257  ft3/day  for  this  one  square  foot  surface  area.  If  the  area  of 
open  water  is  10,000  ft2,  the  amount  of  potential  ice  discharge  is  25,700  ft3 /day,  a staggering 
quantity. 


The  combination  of  frazil  ice  production  and  shore  ice  initiates  and  completes  the  cover  growth 
lor  most  river  systems  in  northern  areas.  Anchor  ice  was  observed  on  the  river  bottom  between 
Sections  5 and  7 where  low-velocity  flows  were  present.  While  one  of  the  authors  was  visiting  this 
site,  several  passing  floes  were  observed  with  debris  attached  to  them,  indicating  that  they  had  sur- 
faced from  the  bed  at  an  upstream  source. 

The  formation  of  shore  ice  and  production  of  frazil  began  in  mid-December  1974  on  the  Con- 
necticut River  and  its  tributaries.  A solid  cover  did  not  remain  in  the  study  reach  until  the  last 
week  of  January  1975,  although  on  several  visits  to  the  Windsor  area,  fragmented  floes  had  arched 
across  the  opening  between  bridge  piers,  establishing  the  formation  of  a temporary  solid  ice  cover 
extending  upstream  to  Sections  5 and  6.  These  ice  covers  were  destroyed  either  by  warm  weather, 
or  by  the  increased  hydrodynamic  forces  caused  by  the  rise  in  river  stage  because  of  the  release  of 
water  for  hydroelectric  generation.  Figure  1 1 gives  an  oblique  view  of  the  arching  of  fragmented 
floes  at  the  bridge  piers. 

The  greatest  potential  for  ice  growth  occurs  over  the  weekends  for  the  Connecticut  River.  Nor- 
mally, no  power  generation  occurs  and  a minimal  flow  is  maintained.  This  constant  water  level 
allows  a uniform  ice  cover  to  form,  but  ice  growth  depends  upon  prevailing  weather. 
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Figure  11.  Oblique  view  of  arching  of  ice  upstream  of  the  bridge  piers,  24  / anuary  1975. 


L 

r 

gss 

r 5 

^ . 

Jf. / . '■  .• 

\ 

S3 

sshi 

When  power  generation  begins  on  Monday,  the  river  rises  and  the  cover  breaks  up  into  a wide 
tange  of  floe  sues.  1 hese  floes  are  then  transported  downstream  where  arching  is  initiated  at  the 
bridge  piers  and  constricted  areas.  The  ice  cover  for  the  1974-75  season  was  observed  to  form  at 
the  bridge  piers  and  propagate  upstream  to  a maximum  point  near  Station  18.  The  thickness  of  the 
ice  was  not  measured,  as  it  was  not  a uniform  cover  of  single  thickness,  but  was  composed  of  frozen 
multiple  thick  floes  of  various  dimensions  and  in  various  orientations. 

(he  first  surveillance  flight  was  flown  on  24  January  1975.  A solid  cover  had  begun  at  the  two 
bridges  in  Windsor  and  was  continuous  upstream  to  about  Station  5.  Upstream  from  Station  6,  a 
fragmented  cover  was  forming  from  floes  that  had  been  broken  as  a result  of  the  power  generation 
that  morning.  The  floes  decreased  in  size  as  they  traveled  downstream,  f igure  1 2 shows  four  se- 
quential oblique  views,  beginning  at  the  covered  bridge  in  Windsor  and  continuing  upstream  to  about 
Section  7.  The  width  of  the  river  in  these  views  is  roughly  500-600  ft.  Frazil  development  was  also 
seen  in  the  reaches  below  Sumner  Falls. 

Ihc  second  and  third  reconnaissance  flights  were  flown  on  7 and  28  February  1975,  respectively. 
On  7 February  solid  ice  extended  to  Section  13  with  a fragmented  cover  to  Section  18.  The  third 
photographic  flight  was  completed  only  two  days  prior  to  the  complete  breakup  of  the  Connecticut 
River  in  the  Windsor  area.  Significant  melt  holes  were  developing  in  the  high  velocity  regions  of  the 
channel,  and  cracking  and  fractured  icc  was  also  visible. 

f igurcs  13  and  14  give  oblique  views  of  the  Windsor  area  on  7 and  28  February  1975,  respectively. 
Effluent  from  fhe  sewage  treatment  plant,  as  well  as  a thermal  effluent,  can  be  easily  seen  just  above 
the  covered  bridge.  Figure  1 5 gives  a view  near  Section  7;  the  melt  hole  to  the  right  is  the  result  of 
both  inflow  Irom  Hubbard  Brook  and  a high  velocity  flow  region  in  the  cross  section.  Another  melt 
region  can  also  be  seen  downstream  along  the  bank  of  the  Vermont  (right)  side,  resulting  from  spring 
runoff  in  the  Runnemcde  Brook. 


N 


19 


a.  Two  views  of  fragmented  ice  upstream  of  the  piers  at  the  coveted  bridge. 


Figure  12.  Four  oblique  views,  beginning  at  covered  bridge  in  Windsor,  and  continuing 
to  about  Section  7,  24  /anuary  1975. 


b.  Upstream  views  of  fragmented  ice  between  Sections  6 and  7. 


Figure  1 2 (cont  d). 


Ice  thickness  profiles  [based  on  an  often-used  empirical  relationship  referred  to  in  a monograph 
by  Michel  (1971)]  were  calculated  using  the  number  of  degree  days  of  frost.  The  formula  is 

V = <x(S)'A 
where  17  = ice  (in.) 

a = coefficient  derived  from  past  experience  (0.7) 

5 = degree  days  of  frost  (°F-days). 

Figure  16  presents  the  potential  static  ice  growth  for  the  last  three  years  based  on  mean  daily  air 
temperature  at  Hanover,  New  Hampshire,  assuming  the  ice  growth  began  on  1 December.  During 
the  winters  of  1972-1973  and  1973-1974,  ice  jams  occurred  in  the  study  reach  between  Sections 
4 and  7.  It  was  suspec^d  that  the  winter  was  much  milder  during  the  1974-1975  season,  but  the 
potential  ice  thickness  growth  appeared  not  to  be  significantly  below  that  of  the  other  two  years 
on  1 March.  The  weekend  air  temperatures  for  these  three  years  were  also  investigated,  and  it  ap- 
pears that  the  complete  freezeup  of  the  river  might  have  occurred  before  mid-January  during  1 972- 
1973  and  1973-1974,  but  no  records  were  available.  The  ice  growth  was  also  calculated,  using 
1 January  1975  as  the  date  of  first  ice  cover  for  the  1974-1975  season,  and  the  results  showed  that 
it  was  below  normal  (see  Fig.  1 6),  but  again  a direct  comparison  cannot  be  made  with  the  other  two 
years. 

Figure  1 7a  shows  the  shearing  and  failure  planes  near  the  shoreline.  The  piling  of  fragmented 
ice  (Fig.  17b)  results  from  the  daily  fluctuation  of  the  river  stage  (3-5  ft).  This  continual  process 
of  shearing  and  piling  is  one  of  the  major  contributors  to  the  formation  of  ice  jams  in  this  reach  of 
the  river.  As  a result  of  the  overnight  freezing,  breaking  of  the  ice  cover,  and  transport  downstream, 
a greater  volume  of  ice  was  produced  and  accumulated  in  the  channel  than  would  result  from  a 
uniform  cover,  if  it  could  be  formed. 

The  sites  of  previous  ice  jams  in  the  Windsor  area  have  ranged  between  the  railroad  bridge  and 
Section  7.  The  shallow  depths,  increased  river  width,  and  sandbars  from  Section  5 to  Section  7 A 
indicate  areas  of  low  velocity.  The  shallow  depth  allows  potential  grounding  of  ice  tloes  and  result- 
ing piling.  The  effect  of  the  backwater  from  the  downstream  Brattlcboro  hydroelectric  dam  on  the 
flow  conditions  in  this  area  has  not  been  established.  The  influence  of  the  solid  ice  cover  in  the 
backwater  could  prohibit  transport  of  ice  from  the  Windsor  reach. 

Figure  1 8 gives  two  views  of  the  study  reach  taken  in  1973.  The  remains  of  the  ice  collar  after 
the  ice  had  receded  about  8 ft  are  still  evident  on  the  middle  bridge  pier.  The  massive  piling  of  ice 
floes  indicates  the  potential  thickness  of  the  jams  and  the  resulting  water  levels  that  can  be  expected 
with  flows  well  below  ice-free  flood  discharges. 


CONCLUSIONS 

The  reasons  that  a jam  did  not  occur  in  the  Windsor  area  during  the  1974-1975  winter  arc  prob- 
ably many:  a solid  ice  cover  did  not  form  until  late  in  January;  the  cover  thickness  appeared  to  be 
less  than  usual;a  mild  breakup  period  occurred  with  no  rain,  leaving  a majority  of  the  ice  to  melt 
in  place;  and  the  White  River  contributed  little  ice. 

A hydraulic  analysis  of  the  flow  conditions  through  the  reach  yielded  certain  criteria  for  identify- 
ing ice  jam  sites.  When  the  surface  conveyance  of  ice  is  significantly  reduced,  the  individual  ice 
floes  tend  to  clog  and  potential  jamming  exists  at  these  sites.  The  results  presented  in  figure  1 1 
reveal  that  the  surface  conveyance  of  ice  improves  with  increased  discharge,  and  that  ice  jams 
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resulting  from  clogging  of  individual  floes  have  a greater  potential  at  low  flows  than  at  high  flows. 

This  simple  ice  discharge  routing  scheme  is  based  on  continuity  of  ice  transport  between  sections. 

No  account  has  been  made  for  the  storage  of  floes,  progression  upstream  of  a cover  if  a jam  is  en- 
countered, or  thickness  of  an  ice  jam  once  formed.  However,  Uzuner  and  Kennedy  (1974)  developed 
a model  for  ice  jam  evolution  to  calculate  the  ice  jam  thickness  and  water  level  from  river  and  ice  char- 
acteristics. Application  of  their  model  should  be  considered  in  future  studies. 

The  Froude  number  analysis  at  each  section  allows  the  determination  of  whether  individual  ice 
blocks  (of  a given  thickness)  will  be  stable  or  will  underturn.  If  the  individual  floes  are  hydrody- 
namically  stable  and  have  sufficient  dimension  and  internal  strength,  arching  of  floes  can  be  ex- 
pected, and  the  cover  can  progress  upstream.  If  the  blocks  are  unstable,  they  will  underturn  and  be 
transported  beneath  the  cover,  with  arching  possible  but  dependent  upon  the  Froude  number.  The 
downstream  distance  and  necessary  hydraulic  conditions  for  block  transport  under  the  ice  cover  have 
yet  to  be  identified,  but  model  studies  are  presently  underway  on  these  topics. 

Ice  jams  on  the  Connecticut  River  in  the  Windsor  area  probably  form  as  a result  of  many  adverse 
conditions,  such  as:  full-depth  constrictions  which  decrease  the  surface  conveyance  of  ice  transport 
and  enhance  the  ice-arching  conditions  (bridge  piers);  potential  blockage  of  ice  transport  near  Sec- 
tion 2 due  to  a solid  cover  at  the  upstream  end  of  the  backwater  from  the  hydroelectric  dam  at 
Brattleboro,  Vermont;  deep  pools  just  upstream  of  the  bridge  piers;  and  shallow  depths  near  Sec- 
tions 5 and  6,  permitting  the  grounding  of  ice  at  low  flows  and  the  increase  of  fragmented  ice  floes  to 
the  area  because  of  the  power  generation  schedule  at  Wilder  Dam. 

The  ice  discharge  parameter  xc  indicates  peak  values  at  the  two  bridges  at  low  flows  (C?  < 5000 
ft3 /sec).  The  flow  Froude  numbers  are  less  than  critical  Froude  numbers  at  these  sections  and  in 
the  immediate  reaches  upstream  of  the  two  bridges  (Sections  4-7a),  enhancing  the  arching  capabili- 
ties of  the  fragmented  ice.  Grounding  of  ice  floes  is  indicated  by  the  critical  Froude  number  being 
equal  to  zero  at  Section  6. 

The  role  of  aerial  photography  in  ice  jam  conditions  cannot  be  overstressed.  Pre-jam  photography 
helped  in  identifying  the  channel  obstructions,  aided  in  surveying  the  cross  sections,  and  provided 
general  information  on  the  channel  characteristics.  The  winter  photography  helped  to  identify  the 
formation  stage  of  ice  cover  growth  and  to  document  the  breakup  of  the  ice. 
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APPENDIX:  LOCATIONS  AND  PLOTS  OF 
CROSS  SECTIONS  FROM  SUMNER  FALLS 
TO  CHASE  ISLAND,  CONNECTICUT  RIVER 
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